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ABSTRACT 

After quasars ionize intergalactic He II at z ~ 3, a large radiation field builds up above the He II 
ionization edge. Unlike the background responsible for H I ionizations, this field should be highly 
variable, thanks to the scarcity of bright quasars and the relatively short attenuation lengths (~ 
50 Mpc) of these high-energy photons. Recent observations of the He II and H I Lyq forests show 
that this background does indeed vary strongly, with substantial fluctuations on scales as small as 
~ 2 Mpc. Here we show that such spatial fluctuation scales are naturally expected in any model 
in which the sources are as rare as bright quasars, so long as the attenuation length is relatively 
small. The correlation length itself is comparable to the attenuation length (> 10 Mpc) for the most 
plausible physical scenarios, but we find order-of-magnitude fluctuations on all scales smaller than 
~ 6 Mpc. Moreover, aliasing along the one-dimensional skewers probed by the He II and H I Lyq 
forests exaggerates these variations, so that order-of-magnitude fluctuations should be observed on 
all scales smaller than ~ 20 Mpc. Complex radiative transfer is therefore not required to explain the 
observed fluctuations, at least at the level of current data. 
Subject headings: cosmology: theory - intergalactic medium - diffuse radiation 



1. INTRODUCTION 



To the vast majority of baryonic matter in the Uni- 
verse, the most important radiation field is the meta- 
galactic ionizing background, and a great deal of ef- 
fort has gone into understanding its properties (e.g. 



Rauch et al l 119 97; McDonald ct al. 2001; Tvtler et al 
20041 iBolton et all 120051: iFan et all 120061: iBecker et al 



2007ft . Most work has focused on measuring the mean 



amplitude of this background, which is indeed the only 
interesting aspect if the background is spatially uniform. 
This appears to be an excellent assumption for H I- 
ionizing photons at z < 4, where spatial fluc tuations are 
probably only a few percent of the mean pud 119923: 
Fardal fc Shull 1993: Mciksin & White 200l lCroftJl2'004 : 



Bolton et all 120051 : iFurlanettol 120081 ). although it breaks 
down nea r and during the rei o nization epoch at 
z > 6 dMeiksin fc White! 12001 IBolton fe HaehneE 
20071 iM esingcr & Diikstra 2008; IChoudhurv et alJl2008t 



Furlanetto fc Mesingerll2008l ). 



Spatial fluctuations are much more important for pho- 
tons above the ionization edge of He II, which can be 
studied with the He II Lyq forest at z ~ 3 (e.g., 
lMiralda-Escudelll993t Dakobsen et all 11994) . A partic- 
ularly powerful approach is to compare the He II and 
H I Lyq forests along a given line of sight. The hard- 
ness ratio rj — Nn e u/N}n parameterizes the strengths 
of each individual absorber; it depends on the ratio of 
the lo cal ionization rates for hydrogen Thi and helium 
ThoII (jMiralda-Escudd fl993h . which are in turn inte- 
grated measures of the ionizing background. Several re- 
cent studies suggest that the He Il-ionizing background 
fluctuates strongly at z ~ 2.6, with nearly order-of- 
magnitude spatial v ariations on scales spanning a few to 
a few tens of Mpc (IShull et al.ll2004 IZheng et ail 12004 
iFechner et al.ll2~006t iFechner fc Reimersll2007al) and some 
evide nce for hardening to higher rcdshifts (Hea p et al.l 
l2000h . 

Large fluctuations are naturally expected for two rea- 



sons. First, only quasars are able to produce photons 
with energies above the 54.4 eV ionization threshold of 
He II, and these sources are relatively rare - implying 
a str ongly fluctuating background even after reioniza- 
tion (|Fardal et all 119981 IBolton et aD 120061 : iFurlanettol 
2008). Direct evidence for these source- induced vari- 
ations has been seen in the "transverse proximity ef- 
fect" of the hardness ratio through comparisons of the 
H I an d He II Lyq forests with surveys for nearby 
quasars dJakobsenet al1l2003l : IWorseck fc Wisotzkl2006t 
iWorseck et all 12007ft . This is exaggerated by the strong 
attenuation from residual intergalactic He II: the mean 
free path of these photons (at > 30 Mpc) is about an 
order of magnitude smaller than that of H I-ionizing 
photons (lHaardt fc Madaul[ l996: Fa ucher-Giguere et al.l 
120081: iFurlanetto fc Ohll2008ft . so a typical point in the 
intergalactic medium (IGM) sees only a few sources. Sec- 
ond, radiative transfer throug h the clumpy IGM can 
induce additional fluctu ations (Masclli fc Ferraral 120051 : 
iTittlev fc Meiksidl2l¥)7ft . 

Naively, one might expect random variations in the 
local quasar population to imprint a characteristic 
scale comparable to the attenuation length, attributing 
smaller scale variations to radiative transfer effects. Here 
we will quantitatively examine the spatial scales of cor- 
relations induced in the He II ionizing background by 
stochastic variations in quasar source counts. We will 
show that the 1 /r 2 intensity profiles around the discrete 
sources create strong small-scale fluctuations and suggest 
that the variations observed in the data may not require 
complex radiative transfer. 

In our numerical calculations, we assume a cosmol- 
ogy with n m = o.26, n A = 0.74, n b = 0.044, h = 

100/ikms- 1 Mpc" 1 (with h = 0.74), n = 0.95, and 
us = 0-8, consistent with the most recent m easurements 
(|Dunklev et alj|2008t iKomatsu et alj|2008ft . Unless oth- 
erwise specified, we use comoving units for all distances. 
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2. THE CORRELATION FUNCTION OF THE IONIZING 
INTENSITY 

We define the correlation function £j(r) of the ampli- 
tude J of the ionizing background (for concreteness, at 
the ionization edge of He II) via 



(•/(ri)J(r 2 )) 



(J> 2 [1 



Ur)}, 



(1) 



where ri : 2 label two points in the IGM, r = |ri — rgl, (..) 
denotes a spatial average, and we have assumed isotropy. 
For Poisson distributed sources and a fixed photon at- 
tenuation length ro, the probability distribution func- 
tion of J(r), as well as the joint distribution function at 
two spat ial points, c an be derived either via Markoff's 
method (IZuol l!992allbT) or the m etho d of charac t eristic 
functions ([Meiksin fc Whltell2003fl ; see lFurlanettol (|2008j ) 
for an a pplication to the He II-ionizing background. IZuol 
(1992b]) shows how to compute the correlation function 
from these distributions: 



(2) 



where Nq = (47r/3)rtQrQ is the mean number of quasars 
inside an attenuation volume, uq is the total number 
density of quasars, and 



dx- 



sinh ; 



exp 



1 



1 - e- 



(3) 



In equation ^ , L denotes the quasar luminosity, and its 
moments are computed by integrating over the source 
luminosity function <&(£). Here we have made several 
simplifying assumptions: (i) quasars counts a re Poisson- 
distributed around the mean value Nq (see IFurlanettol 
2008 for a discussion of deterministic clustering); (ii) 
sources are visible t o infinite di stance, attenuated by a 
factor t = r/r (see lZuolll992bl for a discussion of fluc- 
tuations in attenuation); and (hi) the effective volume 
is Euclidean (certainly reasonable for ro < 300 Mpc, as 
used he re). This corresp onds to the post-reionization 
limit in IFurlanettol (|2008l ). 

The function Ij(u) is shown in IZuol (|1992bf ). It lim- 
its to 7r 2 /4 as u — ► 0, so £j cx r _1 for r -C ro. For 
r ^> r , it is exponentially suppressed. Thus the shape 
of £,/ is determined entirely by our underlying assump- 
tions (and in particular the e~( r / ro ^ /r 2 flux profile of ran- 
domly distributed sources) ; the luminosities and number 
densities of the sources only affect the amplitude of the 
correlation function (although real world complexities, 
and especially radiative transfer, may modify this simple 
expectation). 

However, to determine the amplitude we need (£ 2 ) and 
(L), and thus &(L). This is relatively well-determined 
at the redshifts of interest (2 < z < 3), but unfor- 
tunately our results are extremely sensitive to it. The 
quasar luminosity function is typically parameterized as 
a double powe r law (e.g., iPef fl995t IBovle et al. 2000; 
Hop kins et ai]|2007| ). 



$(L) 



(4) 



where normalizes the density, is the characteristic 
break luminosity, a ~ —1.5, and ~ — 3; all of these 
constants are redshift-dependent. At L 3> L+, $ oc L 13 ] 
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Fig. 1. — Correlation function of the ionizing intensity, £j(r). 
The thick curves use our fiducial quasar luminosity function at 
z = 2.5 and take ro = 35, 55, 110, 350, and 500 Mpc, from top to 
bottom. The thin solid curves show the plausible range of uncer- 
tainty from <t>(-L). 



thus (i 2 ) diverges if j3 > —3 and converges only slowly 
if/3<-3. 

A finite maximum quasar luminosity L max some- 
what alleviates this problem. As a reasonable guess, 
we use the Eddington luminosity for a 10 11 M Q black 
hole (corresponding to a bolometric luminosity 3.3 x 
10 15 L f^), converted t o the i?-band using the correc- 
tion of Hopkins et al. ( 2007) . This is comparable to the 
brightest observed quasars, or ~ 200L* for the redshifts 
of interest. Doubling L max increases the amplitude of £j 
by - 5%} 

Unfortunately, the bright-end slope still introduces 
substantia l uncertainty. To ga uge its impact, consider 
the recent iHopkins et al.1 ()2007t ) luminosity function es- 
timated from a wide range of observational samples. The 
best-fit bolometric luminosity function has « —3.3 at 
z = 2.5. On the other hand, if the fitting parameters are 
forced to vary smoothly with redshift, (3 ss — 3.1. If a 
(non-linear) correction is then applied to reconstruct the 
.B-band values, the luminosity function flattens farther to 
« —2.9. Finally, a naive correction from B-band lumi- 
nosity to the helium ionization edge using the observed 
range of quasar UV spectral indices (jZheng et al.l 119981 
iTelfer et al l l200l iScott et all l200l yields « -2.8 
( Furlanettc 2008f ). although this procedure may overesti- 
mate the flattening if the B-band correction is correlated 
with the far-UV scatter. 

For concreteness, we use the IHopkins et ail (|2007D B- 
band function that varies smoothly with redshift as our 
fiducial model (0 w —2.9) but show some results for other 
plausible values of 0. Note that we only require the total 
number counts N and the shape of the luminosity func- 
tion, so we do not need to correct $(£) to the far-UV 
explicitly. When varying 0, we hold uq constant. 



1 Note that No and the moments of L are also sensitive to the 
minimum quasar luminosity, but this dependence nearly cancels 
from £j. 
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Figure [JJ shows the resulting correlation functions 
at z = 2.5. From top to bottom, the thick curves 
take r = 35, 55, 110, 350, and 500 Mpc. The first 
two are plausible values f or He II-ionizin g photons 
(see iFurlanetto k. Ohl [20081 : IFurlanettol [2008): the oth- 
ers span the range expected for H I-ionizing photons 
(|Madau et al.lll999l : [Faucher-Giguere et al.ll2008fl . 

With Q(L) held fixed between these curves, the am- 
plitude at r -C ro is oc r^/No oc r^ . The shapes, 
when normalized to ro, are identical: £j ~ r _1 at small 
separations followed by an exponential suppression at 
r ~ ro. Measuring this cutoff would therefore be an 
effective measure of the attenuation length. 

However, ro is not the typical scale at which points 
are highly correlated (£j ~ 1). That is set by the num- 
ber density of sources and decreases with ro. We find 
the correlation lengths to be (32, 26, 12, 1.6, 0.8) Mpc for 
r = (35, 55, 110, 350, 500) Mpc, respectively. Only when 
ro ~ 30 Mpc is the correlation length larger than the at- 
tenuation length. 

Instead, the correlation length is closer to (though not 
exactly equal to) the typical radius at which a quasar 
dominates the local photoionization rate, the so-called 
"proximity zone." For He II-ionizing photons, this is 
-Rprox ~ 16 Mpc at z = 2.5 if L B = 10 12 L (near 
Li,) and ThcII ~ 5 x 10~ 15 s _1 (a reasonable value for 
the ionizing backgr ound at z ~ 2.5; Sok asian et alj|2002t 
iBolton et al.l [2006) : -R pr0 x is at least an order of mag- 
nitude smaller for H I-ionizing photons. Strong correla- 
tions are confined to the proximity zones around quasars; 
they are extremely important for the He II photons 
simply because that background flu ctuates so strongly 
([Bolton et alJl2006t lFurlanettdl2008f ). 

The thin solid curves in Figure [1] assume ro = 35 Mpc 
but take different values for (3 in &(L): the upper curve 
has /3 « —2.8 (corresponding to a naive inclusion of the 
scatter in quasar UV spectral; IFurlanettol l2008ft . while 
the lower curve follows the shape of t he z — 2.5 (only) 
bolometric fit in lHopkins et all (|2007[ ). with /3 « -3.3. 
The total uncertainty in the amplitude of £j therefore 
spans at least a factor of three, comparable to the differ- 
ences from the plausible range of ro . 

Note that we have assumed here that quasars pro- 
duce the entire ionizing background. While most likely 
true for He II-ionizing photons (though see the discus- 
sion of recombination radiation in §4.2p , galaxies provide 
much (if not most) of the H I-ionizing background (e.g., 
iFaucher-Giguere et alj |2008| ) . Ga laxies are so com mon 
that their fluctuations will be tiny (jFu rlanetto 2008]) , and 
the additional background can be considered uniform. In 
that case, £j is suppressed by a factor (JQso/^tot) 2 , and 
the correlation length will decrease as well. 



3. THE POWER SPECTRUM OF THE IONIZING 
INTENSITY 

A cleaner tool for separating the scale-dependence of 
the correlations is the (three-dimensional) Fourier trans- 
form of the correlation function, the power spectrum 




k (comoving Mpc -1 ) 

Fig. 2. — Three-dimensional power spectrum of the ion- 
izing intensity, normalized to dimensionless form Aj(fc) = 
[(fc 3 /27r 2 )Pj(fc)] 1 / 2 . All thick curves assume z = 2.5 and take 
ro = 35, 55, 110, 350, and 500 Mpc, from top to bottom (line styles 
are identical to Fig. [TJ; the thin solid curves show the plausible 
range of uncertainty from &(L). 



P.j{k): 2 

Pj(k)=^fdrUr) r -^^-. (5) 

We show the dimensionless form of this quantity, 
Aj(k) = [(fc 3 /2 7 r 2 )Pj(fc)] 1 / 2 , in Figured this is roughly 
equal to the fractional rms scatter when the ionizing 
background is smoothed on a scale i = 2n/k. The line 
styles are identical to those in Figure [TJ with ro increas- 
ing from top to bottom. 

As one should expect, Aj has a break at fco = 2-k/tq. 
At k > fco, A j cx /c 1 / 2 , which follows from the 1/r de- 
pendence of £j. For k < ko, Aj cx fc 3//2 , thanks to the 
exponential suppression in £j. 

Unfortunately, measuring this power spectrum, and 
hence identifying this feature, requires three-dimensional 
data. For studying the ionizing background, we are in- 
stead limited to one-dimensional skewers along the rare 
Lya forest lines of sight to UV- bright quasars. Un- 
der these conditions, we can only measure the one- 
dimensional power spectrum, 

PjMH) = J k £kPj(k), (6) 

where the integration over smaller physical scales ac- 
counts for the possibility of aliasing. 

This immediately presents a problem: at k > fco, Pj oc 
fc~ 2 , so the integral is logarithmically divergent. Clearly 
the observed power at every scale will be dominated by 
the smallest scales for which the power is non-zero. In 

2 Actually, a s impler way to compute the power spectrum is with 
the halo model (Cooray & Shcth 2002); our physical model of un- 
corrected sources with constant attenuation is equivalent to ignor- 
ing the two-halo term and inse rting a "halo" profile cx e~ ( r / r o) / r 2 
llScherrer k. Bcrtschingcr 19911'). 
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Fig. 3. — Fractional rms scatter in the radiation field along 
one-dimensional skewers, as a function of smoothing scale (or 
wavcnumber k). All thick curves assume z = 2.5 and take 
ro = 35, 55, 110, 350, and 500 Mpc, from top to bottom (line styles 
are identical to Fig. [TJ; the thin solid curves show the plausible 
range of uncertainty from <&(!/). 



reality, this limit will be set by Jeans smoothing and 
thermal broadening in the IGM, which together smear 
out power on smaller scales. For concreteness, we set 
fcmax = 10 Mpc^ 1 in our calculations, comparable to the 
Jeans scale for fully- ionized mean density gas at z = 2.5. 

Figure [3] shows the one-dimensional power spectrum 
of the ionizing background, again presented in non- 
dimensional form such that the amplitude is the rms scat- 
ter expected when data are smoothed on a given physi- 
cal scale. The curves are identical to those in Figure O 
The exponential suppression at k — 10 Mpc -1 reflects 
the assumed Jeans scale. Because of aliasing, the one- 
dimensional variance is always dominated by the smallest 
scales, and there is only a slight break at ro in Pj t iT>(k); 
it will unfortunately be very difficult to extract this pa- 
rameter from one-dimensional Lya forest data. 

However, the amplitude of Pj^u does depend strongly 
on ro and offers another route to estimate it, albeit with 
relatively large uncertainties. Unfortunately, the plausi- 
ble range of <£>(£) (illustrated by the thin solid curves), 
introduces about a factor of two uncertainty in the over- 
all amplitude and will make it difficult to measure ro 
robustly. 

We also find that the fluctuation amplitudes are quite 
large: for ro = 35 Mpc, the fractional variations are of 
order unity even across hundreds of Mpc, and an order 
of magnitude at scales ~ 20 Mpc. In contrast, the varia- 
tions in the three-dimensional field only reach this level 
at ~ 6 Mpc. Even when ro — 500 Mpc, order unity fluc- 
tuations appear at ~ 10 Mpc. 

4. DISCUSSION 

We have shown that the scarcity of quasars induces 
strong correlations in the radiation background, even if 
those sources are themselves randomly distributed. This 
is especially important for He II-ionizing photons, whose 
attenuation lengths are < 50 Mpc after He II reioniza- 



tion, and so have expected correlation lengths ^25 Mpc. 
The variations are especially large when viewed along 
one-dimensional skewers (as provided by the Lya forest): 
in those circumstances aliasing induces order-unity fluc- 
tuations on scales > 100 Mpc, and order-of-magnitude 
fluctuations on scales < 20 Mpc. 

4.1. Comparison to Observations and Simulations 

The best way to observe these correlations is by com- 
paring the He II and H I Lya forests. The relative abun- 
dance of these species depends primarily on the relative 
radiation backgrounds above their ionization thresholds, 
so the hardness parameter rj = -/VhcIi/^Vhi provides a 
measure of the local str e ngth of the He II i o nization rate 
(iMiralda -Escudd ll993UShull et all 12004 [Zheng et all 
120041 iFechner et al.ll200"6tlFechner k. Reimeril2007al ). In 
the limit of a uniform Thi, and if the line comparison can 
be made cleanly, rj ex 

l/r He ii- 3 In that case, we can com- 
pute £j directly from the observations because £iy, ? = £j 
(IZuolll992bD . 

Although this correlation function has not yet been 
measured precisely, the data clearly show l arge 77 vari- 
ations on rather small scales. For example, [Shull et al.l 
(|2004f ) detected substantial r\ fluctuations over ~ 2 Mpc 
segments toward HE 2347-4342. Bu t their quantitative 
level re mains unclear. For example, Fec hner fc Reimersl 
(|2007al) showed that the highest signal-to-noise portions 
of that spectrum, and the entire line of sight to HS 
1700+6416 (also high signal-to- noise) , are dominated by 
smoother variations, with ~ 33% of the IGM varying 
on scales < 6 Mpc and > 50% varying only on scales 
> 14 Mpc. 

Contrary to previous claims, these qualitative measure- 
ments are consistent with the large variations on mod- 
erate and small scales expected from a simple model 
of discrete quasar sources, even though the attenuation 
length is much larger. In particular, we emphasize that 
aliasing amplifies the apparent fluctuations along Lya 
forest skewers, so other processes - such as complex 
radiative transfer - may not be required to reproduce 
the observations ( although they may still be im p or t ant ; 
iMaselli fc Ferrarall2005t iTittrev k. Meiksinll2007h . 

Quasar-induced fluctuations have been directly de- 
tected through the "transverse p roximity effect" in ri 
along two separate lines of sight (I Jakobsen et alJ 120031 : 
iWorseck fc Wisotzkil [20061 iWorseck et al.ll2007ft : the ra- 
diation background hardens near several quasars close to 
the line of sight. However, even within the regions of low 
77, there are still large fluctuations, especially at the lower 
rcdshifts (where ro is larger), and there are other regions 
with low 77 but no sources within ~ 30 Mpc. Whether 
these differences are due to aliasing, anisotropic emission, 
or radiative transfer remains to be seen. 

Although detailed numerical simulations have not yet 
addressed thi s question, su p port f or our conclusion was 
provided by iBolton et al.l (|2006f ). who found strong 

3 We should note that it is not trivial to comp are lines in this way, 
especi ally if they are thermally broadened ||R3chncr &; Reimersl 
2007b), and in some cases the "pixel optical depth method" may 
be sup erior. Fortunately, the two methods appear to agree in most 
cases (Fechner & Rcimers 2007a). Direct comparisons to numeri- 
cal simulations is probably the safest route, but it has not yet been 
possible (largely because of the difficulty of properly simu lating 
Hell reionization; IPaschos et a l. 2007; McQuinn ct al. 200|). 
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small-scale variations in simulated data using a model 
of discrete quasar sources without appealing to true ra- 
diative transfer (see their Fig. 5). The next step is to 
compute the power spectrum of the data in order to 
quantitatively compare to models like ours. Whether 
that is possible with the sparse and noisy data currently 
available remains to be seen, but many more lines of 
sight have recently become available (|Zheng et al.ll2008t 
ISvphers et al.ll2008l ). so these tests should sharpen in the 
near future. 

4.2. Caveats to the Model 

However, at the moment our model itself is also too 
crude to make detailed comparisons with the data, and 
improvements are clearly needed on the theoretical side. 
In particular, several effects will conspire to decrease the 
observed fluctuations relative to our calculations. First, 
we have included the entire intensity profile around each 
source, including the (divergent) r — > limit. Of course, 
structure within the quasars host halo will actually mod- 
ify the small-scale profile (and it is not measured by the 
Lya forest anyway), decreasing some of the high-fc power. 

More importantly, we have ignored the limited dy- 
namic range available in measurements of r/: in reality, 
we cannot measure either extremely large values of r\ (be- 
cause of saturated absorption) or extremely small values 
(where the transmission is essentially unity). The lat- 
ter limit is especially important, because it will damp 
the observed power within the highly- transpa rent central 
proximity zone of each source. IZuol (|1992bl ) considered 
some of the complexities of the nonlinear transforma- 
ti on, but this ef f ect is best studied with mock spectra as 
in lBolton et~aT1 (|2006D . 

Third, we have only considered photons with a sin- 
gle, well-specified attenuation length (and chosen ro for 
photons just above the He II ionization edge). However, 
quasars ha ve hard spectra (ty pically with L v cx in 
the far-UV; iTelfer et al.ll2002f ) and so produce a substan- 
tial number of high-energy photons. Because such pho- 
tons ha ve longer mean fr ee paths, they damp the corre- 
lations (jFurlanettol 2008). However, the photoionization 
cross section falls steeply with frequency (cr v oc so 
this should not qualitatively change our conclusions. 

Another concern is our assumption that the en- 
tire He II-ionizing background comes from quasars. 
While galaxies are un likely to be significant (but see 
iFurlanetto fc Of] [2008). a diffuse background must de- 
velop from the recombination radiation that follows IGM 
ionizations. The importance of this background depends 
on the geometry of the absorbers: provided that most are 
optically thick regions, much of the recombination radi- 
ation will be directed deeper into the (mostly neutral) 
syste m and thus not influence the IGM (jMiralda- Escudel 
l2003h . However, in the opposite case-B limit, ~ 40% of 
the photons are re-emitted. This could therefore decrease 
the amplitude of £j by a factor of order unity. 

One might also wonder whether the assumption of a 
uniform Thi is acceptable for these purposes: inside the 
proximity zone of each quasar, both ThcII and Thi in- 
crease, such that r\ may remain constant. However, we 
have shown that, for the large attenuation lengths rele- 
vant to H I-ionizing photons, the correlations induced by 
discrete quasars only manifest themselves on < 1 Mpc 
scales, and this will be reduced still farther if galaxies 



dominate the H I-ionizing background (as now seems 
likely; e.g.. lFaucher-Giguere et~atll2008h . This is over an 
order of magnitude smaller than the correlation length 
for He II-ionizing photons. In other words, because so 
many more sources contribute to the local H I-ionizing 
field, the H I proximity zone is much smaller than that 
for He II, so it cannot compensate for the fluctuations in 
the higher-energy background. 

On the other hand, we have ignored the deterministic 
bias of quasar hosts and included the stro ng variations 
seen in far-UV quasar spectra only crudely (|Telfer et al.1 
l2002tlScott et al.l2004h . Both of these effects will amplify 
the small-scale correlations, although the degree remains 
uncertain. 

The large fluctuations we have found may seem sur- 
prising in comp a rison to previous studies, such as 
IMeiksin fc White! (|2003t) . They considered a quasar- 
dominated ionizing background at high redshifts [z > 4) , 
where the fluctuations in J are also large. But they 
found that their effects on the H I Lya forest trans- 
mission statistics (such as the flux probability distribu- 
tion function and its power spectrum) remain modest. 
Several aspects of those statistics probably account for 
the difference with our conclusion. First, the Lya opti- 
cal depth, on its own, is most sensitive to IGM density 
fluctuations, which may mask variations in th e radiation 
background (e.g.. IWorseck fc Wisotzkil 12006 .) . Second, 
the complex line structure of the forest (peculiar veloc- 
ities, line broadening, etc.) may conceal some of the 
r fluctuations. Finally, the e~ T suppression decreases 
the available dynamic range in the forest, as described 
above. Direct comparison of the He II and H I forests 
provides a cleaner tool to study these fluctuations prin- 
cipally because it removes the degeneracy with density 
(which affects both lines t he sam e way, at least in prin- 
ciple). IMeiksin fc White! (|2003f ) also assumed a signifi- 
cantly steeper luminosity function (/3 = —3.41) than ours 
(/3 « —2.9), which reduces the fluctuations by a factor 
~ 2. 

4.3. Future Directions 

In addition to the improvements outlined above, spa- 
tial correlations in the high-energy radiation field can 
be used to address two other sets of questions. First, 
we have only calculated £j after He II reionization has 
ended. During reionization, t he ionizing back ground 
fluctuates even more strongly (|Furlanettol |200 8\ with 
fewer discrete sources illuminating any point in the IGM. 
Thus the correlations will be even stronger - and, indeed, 
in many ionized regions there will only be a single active 
quasar, which could allow us to trace out its region of 
influence in detail. 

Second, to the extent that these correlations reflect the 
proximity zones of individual bright quasars, we may be 
able to pick out regions with large ThoU from the ob- 
served spectra and associate them with known quasar 
(or at least galaxy) locations. In the best case, we could 
trace the intensity profile around each nearby source 
and so test for the "transverse proximity effect;" such 
tests have already shown hardening in the local radia- 
tion field near quasars and placed bounds on the life- 
times of UV-bright quasars (jWorseck fc Wisotzkil 120061 : 
IWorseck et al.l I2007D . In the crudest sense, such re- 
gions would suggest target regions for future quasar 



(> 



searches (as in lJakobsen et all |200 3f). These studies re- 
quire high signal-to-noise measurements of the Lya for- 
est (in both He II and H I) and also detailed maps of 
their environments, ideally spanning a couple of attenu- 
ation lengths (or ~ 100 M p c, several times farth er than 
iWorseck fc Wisotzkil [20M iWorseck etHl 12001 . This 
exercise would elucidate, in detail, not only the attenua- 
tion of high-energy photons throughout the IGM but also 
properties of quasar host galaxies, the lifetime and light 
curves of quasar sources (by tracing the intensity profiles 



of each source) , and the geometry of the quasar emission 
(i.e., whether it is isotropic and whether it evolves over 
the lifetime of the quasar). 
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